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ABSTRACT 

The dynamics model and data sources used herein 
to perform air-data reconstruction are discussed, 
and a brief discussion of the Kalman filter is 
Included. The discussion indicates the need for 
adaptive determination of the noise statistics of 
the process. The filter innovations are presented 
as a means of developing the adaptive criterion, 
which is based on the estimation of the true mean 
and covariance of tho filter innovations. A method 
for -he numerical approximation of tho moan and 
covariance innovations is presented. 

The algorithm as developed is applied to air- 
data reconstruction for the space shuttle, and data 
obtained from the third landing are presented. To 
verify the performance of the adaptive algorithm, 
the reconstruction is also performed using a 
constant covariance Kalman filter. The results of 
the reconstructions are compared, and the adaptive 
algorithm exhibits better performance. 


SYMBOLS 

A kinematic acceleration of aircraft 

b white, GausBian distributed unity 

covariance error in measured system 
response 

C system geometry matrix, L x M dimen- 

sioned matrix 

e filter innovations vector, L x 1 

dimensioned 

F state-noise gain vector, M x 1 dimen- 

sioned 


FF discreto-time state-noise covariance, 

M x M dimensioned 

f rolloff frequency of low-pass filter 

used in dynamics equation 

G measurement noise vector, L x 1 

dimensioned 


T 

GG moaouroment noise covariance matrix, 

L x L dimensioned 

k quantity corresponding to kth sample 

L dimensions of measurement vector 

M dimensions of state vector 

n nonstationary disturbance to state 

dynamics 

P filter error covariance matrix, M x M 

dimensioned 

t continuous-time index 

V groundspeed state used in dynamics 

equation 

Vw airspeed state used in dynamics 

equation 

W windepeed state used in dynamics 

equation 

X position state used in dynamics 

equation 

x state vector, M dimensioned 

z measurement vector, L dimensioned 

4 discrete transition matrix, M x m 

dimensioned 

Mathematical operators! 

COV( ) covariance of quantity ( ) 

d/dt( } derivative of ( > with respect to 

time 

E( ) expectation of quantity { ) 

S( ) sample mean of quantity ( ) 

T 

( ) transpose of ( ) 

( ) quantity es'timated by filter 
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Subscripts t 

k value of quantity at sample time k 

k+1/k indicates quantity resulting from 

extrapolation atop of filter 

k+1/k+1 indicates quantity resulting from 

correction step of filter 

INTRODUCTION 

Accurate estimates of the air data of high- 
speed aircraft are essential to successfully per- 
form flight research and developmental tusks • 

Recent expansions of operational flight envelopes 
have necessitated the development of advancod 
reconstruction techniques to obtain these esti- 
mates > Many of these techniques use a Kalman 
filter to merge several data sources to enhance the 
result. Most of these methods assume that random 
disturbances in the aircraft dynamics can be 
accounted for by assuming constant filter statin- 
tics which are assumed to be known a priori . 
Unfortunately, thore is often little physical evi- 
dence by which one can determine the statistics of 
those random disturbances. In many cases the dis- 
turbances are nonstationary) in other words, in 
such cases the random disturbances in the aircraft 
dynamics can only be described by time- varying sta- 
tistics. These cases include such unsteady regimes 
as flight in turbulence, high-angle-of-attack 
flight, and aerostructural buffet. In such cases, 
the statistics of random disturbances are often 
selected purely by guesswork. The so-called "opti- 
mal" results of ouch filters are optimal only in 
the sense that an error criterion is minimised. 

For nonstationary cases, the estimation scheme 
must allow for tfre adaptive estimation of the sta- 
tistics that describe random disturbances in the 
aircraft dynamics as part of the filter loop. This 
paper presents such an algorithm. The algorithm, 
although primarily intended for applications to 
air-data reconstruction, has applications to a 
variety of other fields. Because the space shuttle 
ie known to encounter several or all of the above- 
mentioned unsteady flight conditions during its 
reentry, space shuttle reentry data are used to 
verity the resulting algorithm. 

ESTIMATION CONCEPT 

The estimation problem is essentially one of 
complementary filtering. Data from four indepen- 
dent measurement sources are merged by means of the 
filtering algorithm to give an enhanced result. 

The data sources include high-, medium-, and low- 
frequency data. The high-frequency data are pro- 
vided by a s trap down linear accelerometer package) 
the medium-frequency data are provided by a pair of 
pneumatic hemispherical air-data seneors and C-band 
radar tracking) and the low-frequency data are pro- 
vided by a meteorological analysis of the atmos- 
pheric conditions. The air data are complementary 
to and blended with the meteorological data) the 
acceleration data are complementary to and blended 


with t*>e tracking data. Those are in turn com- 
bined to give enhanced estimates of the aircraft 
poeitlon, groundspeod, windopoed, and airepeed. 

The resulting estimates possess characteristics of 
all four measurement sourcos. A schematic of this 
concept is presented in Fig. 1. Each of the data 
sources is described in detail In the INFORMATION 
SOURCES section. 

INFORMATION BOURCES 

Tracking Data 

t 

Radar tracking data are obtained from a akin 
track using an FPS-16, C-band, hlgh-rango tracking 
system (Ref. 1). Provided are highly accurate, 
medium-frequency measurements of the aircraft's 
range, azimuth, and elevation relative to the radar 
site. The tracking data were recorded at 20 Hz and 
interpolated to 25 Hz. 

Meteorological Data 

Steady-state meteorologically derived wind and 
barometric data are used to provide accurate but 
very low frequency information concerning the dyna- 
mics of the atmosphere along the reentry flight- 
path. The data are obtained by a series of weather 
balloons launched at various times and locations 
along the anticipated flightpath. The raw data 
thus obtained wore corrected for diurnal and spa- 
tial variations (Ref. 2). Examination of time and 
altitude variations in the data gives indications 
of both the steady-state magnitude and the turbul- 
ence in the winds aloft. The motorological data 
were interpolated to 25 Hz using radar position 
data. 

Air-Data Measurements 

The air-data measurements are obtained from 
a pair of side-mounted hemispherical pneumatic 
sensor (Ref. 2). In addition to sensing turbu- 
lence and compression caused by the local flow 
field, the sensors are subject to pneumatic lag. 

These factors are manifested as nonstationary 
disturbances. The hemispherical sensor data are 
recorded at 12.5 Hz and interpolated to 25 Hz. 

Strapdown linear Accelerometer Data 

High-frequency data concerning the aircraft's 
inertial dynamics are provided by an onboard, 
strapdown, linear accelerometer package (Ref. 2). 

The accelerometer package, intended for aerodynamic 
coefficient identification, provides very high fre- 
quency measurements. The package, however, is not * 

inertial quality and is subject to significant bias 
errors. Thus resulting data cannot be integrated 
open loop. The strapdown data wore rotated to 
earth-relative topodetic coordinates before use in 'f 

the filtgr. The direction cosines used in perform- 
ing the rotation were obtained directly from the 
shuttle's Inertial measurement unit (Ref. 3). 

After rotation, the acceleration of gravity as a 
function of' altitude was added to the vertical 
component. The strapdown data are recorded at 
176.0 Hz and decimated to 25 Hz. 
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FILTER EQUATIONS 


Tho discrete-time form of tho Kalman filter 
was ubo<1 an the starting point for developing the 
algorithm. Tho filter an mechanized conniata of 
threo parallel filters, one for each topodotic axis 
component. The topodotic axis system (Fig. 2) is 
defined no as to have its x-axia directed northward 
with reepeot to tho local horizon, the y-axis 
directed eastward, and the z-axia directad toward 
tho center of the earth. Each component filter was 
assumed to have four states i aircraft position, 
groundspeed, windspeod, and true airspeed. The 
dynamics equations, although continuous time, are 
easily discretized by transition matrix integra- 
tion. Tho matrix equation chosen to describe the 
process dynamics is 
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The relationships of the states to the measured 
data are 
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+ G(t) b(t)c 


In Eqs. (la) and (1b), tho parameter f repre- 
sents a first-order low-pass filter rolloff fre- 
quency. This parameter affects only the windspeod 
and airspeed components of the filter. The low- 
pass filter terms are chosen so as to insure that 
the resulting windspeed and airspeed terms will be 
smooth. 

The function n(t) represents a nonstationary 
random disturbance in the proceB3 at time t. It 
is assumed that the model of Eq. (la) adequately 
accounts for all systematic aircraft dynamics! as 
a result, n(t) is assumed to be zero mean. The 
vector F{ t) , a four-element vector, is used to 
describe the relative effect of random disturbances 
occurring at time t on eaoh of the states. Since 
n(t) is assumed to be zero mean, the vector product 
of F with itself is assumed to approximate the 
covariance of tho random disturbances at time t. 
Consequently, it is the vector product of F with 
itself that must be identified by moans of the 
adaptive criterion developed in the ADAPTATION 
CRITERION section. 

The function b(t) represents an error in the 
measurement vector at time b. This quantity is 
assumed to he locally stationary (that is, it has 
statistics that vary slowly with time), and it has 
a nonzero mean. The vector G(t), a three-element 
vector, is used to describe the respective ampli- 
tudes of the errors in each measurement at time t. 
The elements of G(t) are essentially measurement 


weighting factors. As a part of tho adaptivo 
algorithm presented in tho ESTIMATION OF INNO- 
VATIONS MEAN soction, tho bias errors in tho 
measured data are estimated and compensated for. 

For this reason, the vector product of G(t) with 
Itself is assumed to approximate the covariance of 
the measurement orrors at timo t. 

The first olomont of G(t) 1 b used to weight the 
erroru that are expacted to occur in the radar- 
derived position data. Radar position orrors are 
known to become significant for elevation anglos of 
less than 10 s . For this reason the first olomont 
of G(t) is prescribed to dccreaee linearly with 
increasing elevation angle until elevation angle 
reaches 10*. At elevation angles greater than 10", 
the first element of G(t) is hold constant. For 
elevation angles less than 0 s , tho value of tho 
first element of G(t) is fixed at 1. A schematic 
of this weighting scheme is presented in Fig. 3(a). 

The second element of G(t) is used to weight 
the expected errors in the hemispherical sensor- 
derived airspeed data. Hemispherical sensor data 
are known to become poor at high angles of attack. 
Consequently, for angles of attack greater than 
10", the second element of G(t) is proscribed to 
increase linearly with angle of attack. For anglos 
of attack less than 10", the weighting is held 
constant at 1. A schematic of this weighting 
scheme is depicted in Fig. (3b). 

The third element of G(t) is used to weight the 
expected errors in the meteorologically derived 
windspeed data. Because little precise information 
ie available regarding the accuracy of these 
measurments, those data are considered to be 
equally accurate throughout the flight envelopo. 

For this reason, tho third element of G(t) is held 
constant at 1. 

ADAPTATION CRITERION 

As mentioned previously, tho statistics of the 
random disturbances in the process dynamics are 
assumed to vary as a function of time. For this 
reason, wo must develop an adaptive criterion by 
which these statistics can be estimated. Tho adap- 
tive scheme discussed in thiB paper uses infor- 
mation provided by the filter error vector. The 
error vector is defined as the difference between 
the measured system response and the expected 
Bystem response. This vector, 0^+1, is also called 
the filter innovations, which is the name used 
henceforth. The statistics of the filter innova- 
tions tell a great deal about how well the filter 
model is performing. Because of tho assumptions 
UBed in deriving the discrete Kalman filter, the 
innovations should be a purely white, Gausslan- 
distributed, zero-mean sequence. If such is the 
case, then one can be fairly confident that the 
filter model is correct and that the resulting 
estimates are close to optimal (Ref, 4). If this 
is not the case, then the parameters of the filter 
model are in error and information provided by the 
innovations can be used to drive the assumed model 
toward the correct model. This process is depicted 
in Fig. 4. 
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ESTIMATION OP INNOVATIONS MEAN 

The Kalman filter is an unbiased estimator; 
consequently, for unbiased measurements, the moan 
of the innovations should be exactly zero (Ref. 5). 
If, however, there are biases in the measurements, 
the moan of the innovations will equal the expected 
bias orror in tho measurements. Measurement bias 
errors do not usually change rapidly with time, and 
for this reason the true moan can be approximated 
by sampled statistics. The sample mean can be 
allowed to change slowly with time by taking the 
time average over only a fairly local region. Tho 
time-averaged estimate of the bias error s{o k |, can 
bo recursively subtracted from the measurement vec- 
tor to form a transformed measurement vector that 
is nearly unbiased. Using this transformed vector, 
the estimation algorithm proceeds us in the stan- 
dard Kalman filter. 

ESTIMATION OF THE INNOVATIONS COVARIANCE 

The negation of the approximate biases from tho 
measurement vector allows for considerable simpli- 
fication of the adaptive process; the problem redu- 
ces to one of estimating the mean square of the 
innovations. This quantity is estimated in the 
following manner; The actual system response is 
given to be related to the system states by some 
geometry matrix, C, and an error term, b(t); that 
is, 

=k “ c x k + =k b k < 2b > 

Therefore, the oxpectod system response is given by 

z k » Cx k (2b) 

As a result (remembering that biases have been 
removed), the innovations covariance is described 
by 

COV(e k ) “ E(z k - Cx k/k . t )(z k ^ Cx k/k „ 1 ) T (3a) 
Substituting in Eq. (2a) for z gives 

COV(e k ) - E[c{x k - X k/k _,) + G k n k ][c(x k 

" x kA-l) + °k n k] T 

Now, since the measurement noise is assumed to be 
white and all biases have been previously removed, 
the previous expression reduces to 

COV(e k ) - CE(x k - x k / k -i) (x k - x k/k _ 1 ) T C T + Gg£ 

From the definition for the filter error covariance 
matrix, 

p k/k-1 " E ( x k “ X k/K~l) ( x k “ x kA-l) T 

we get 

COv(e k ) - CP^-! C T + GGjJ (3b) 


The quantity COV e k is computed as a part of 
the gain expression in tho standard implementation 
of the discrete Kalman filtor. Conveniently, using 
the previously statod assumptions, tho covariance 
of the innovation is available for use with no 
additional computational expense. 

ADAPTIVE COMPUTATION OF THE 
STATE-NOISE COVARIANCE 

For the standard formulation of the dlBarote 
Kalman filter, the state-noise covariance extrap- 
olation stop is given by 

P k+1A “ * p kA + FF k <*«> 

T 

whore FF k is the state-noise covariance and 4 
is the system transition matrix. Substituting 
Eq. (4a) into Eq. (3b) and rearranging gives 

CFf£c T - COV(c k ) - c(4P k . 1/kM 4 T )c - GG k (4b) 

If tho geometry matrix, C k » is square and non- 
singular, we can t solve for the state-noise covari- 
ance matrix directly. If the number of measure- 
ments does not equal the number of states, we can 
expand the matrix equation to form a sot of L 
scalar equations in M(M - 1) unknowns, whore L is 
the number of measurements and M is the number of 
states* Wo then must assume that the off-diagonal 
elements of the covariance matrix are insignifi- 
cant, and the covariance matrix is essentially 
diagonal. This allows us to solve directly for l> 
of the diagonal elements in the state-noise covari- 
ance matrix. In the case of Eq. (1), we can solve 
directly for the first, second, and fourth diagonal 
elements; the remaining element can be solved for 
by considering that by definition 

Vw - V + W (4c) 

from which we reason that the third diagonal ele- 
ment is simply the cum of the second and fourth 
diagonal elements. Equation (4b) can now be imple- 
mented as a part of the filter loop to give a 
closed-loop estimate of the covariance of random 
disturbances in the aircraft dynamics. 

PRESENTATION OF ADAPTIVE TECHNIQUE 

The adaptive algorithm is presented in schema- 
tic form in Fig. 5. The portion of Fig. 5 that 
lies outside the dashed line represents the procoss 
dynamics; the portion that lies inside the dashed 
line represents the estimation loop. The arrows 
depict the flow of information through the filter. 
The state-noise covariance is computed by means of 
Eq. 4 as a part of the recursion of the filter. 

APPLICATION OF ALGORITHM TO TRAJECTORY 
RECONSTRUCTION PROBLEM 

As a verification of its validity, the adaptive 
algorithm is now applied to the problem of air-data 
reconstruction. Data obtained from the third space 
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shuttle reentry (STS-3) era chosen to illuatrata 
the problem. The STS-3 landing at White Sands, Hew 
Mexico, occurred on a worst-case day. Extremely 
high winds and moderate-to-sovore turbulence wore 
known to exist. Examination of landing-day rawin- 
sonde and Jimsphere balloon data provided by tha 
Air Force Flight Test Center at Edwards APB indi- 
cated that jetatroam velocities were measured in 
excess of 150 knots. Tho Jimsphere balloon was 
observed to have rise-rato oscillations that varied 
from 0 to 10 ft/sec. Those wind conditions 
extended throughout the troposphere and the lower 
stratoBphoro. Because of these conditions, tho 
vertical wind component was subject to greater ran- 
dom variation than has bean experienced during any 
of the other spate shuttle roontrioe. Tho condi- 
tions are suggestive of mountain gravity wave 
activity. Analytical solutions documented in 
Ref. 2 also suggost that this is the case. Such 
conditions are known to produco highly nonstation- 
ary disturbances (Refs. 2 and 6). 

The air-data reconstruction was performed 
using the adaptive algorithm with arbitrary initial 
covariances assumed for the state disturbances. 

Tho measurement error covariances wore chosen based 
upon the prescribed methods mentioned earlier. To 
verify the effectiveness of the adaptive algorithm, 
the air-data reconstruction was also performed 
using s standard implementation of the discrete 
Kalman filter. The standard implementation was 
performed using the same initial covariances as 
abovo; however, in this case both the state distur- 
bance and measurement error covariances were held 
constant. 

Tha results of the adaptive case are presented 
in Fig. 6. Presented are comparisons of the 
measured and reconstructed time histories of air- 
speed {Fig. 6(a)), vertical windspeed (Fig. 6(b)), 
Jnd total horizontal windspeed (Fig. (6c)). Esti- 
mates of these quantities are tho most heavily 
affected by time variations in random atmospheric 
disturbances. Inertial type quantities such as 
position and groundspeed are less hoavily affected. 
The results show no significant discrepancies 
between tho measured and estimated values. Con- 
sidering the fact that mountain wave activity was 
believed to exist, the vertical windspeeds, 
although large, are within believable limits. Ab 
expected, the filtered estimates exhibit higher 
frequencies than do the meteorological winds. 

The results of the nonadaptive estimates are 
presented in Fig. 7, and similar comparisons are 
made. These comparisons show significant discre- 
pancies. The estimates of airspeed and horizontal 
windspeed differ by more than 100 ft/sec. It is 
not believed that the meteorological estimates 
could have been this much in error. The verti- 
cal wind estimate reaches a peak value of nearly 
60 ft/sec downdraft'. Had these conditions actually 
occurred, it is doubtful that the reentering space 
shuttle could have cleared the lofty San Andreas 


Mountains that border tho STS-3 landing site. In 
Fig. 7 it is interesting to note that the behavior 
of the standard filter becomes more roaliauic 
toward the and of the tima histories. Thosu data 
correspond to data that were obtained just before 
landing and several minutes aftor the shuttla had 
dropped below the regions of heavy turbulance and 
mountain wave activity. Under these conditions the 
disturbances dropped to nearly zero and no adaptive 
estimation of the covariances of the disturbances 
wore necessary. 

CONCLUDING REMARKS 

An adaptive algorithm that can bo used to 
estimate the state-noise covariance for certain 
typos of nonstationary processes has boon deve- 
loped. The algorithm, which was developed pri- 
marily for tho purpose of air-data reconstruction, 
accounts for improper knowledge of the state- 
disturbance covariance matrix and, to some oxtont, 
accounts for unknown msasuroment biases. Tho algo- 
rithm is recursive and has the potential for real- 
time implementation. Because the adaptation cri- 
terion wsb formulated in a general sense, it has 
applications to fields other than air-data 
reconstruction. 

An air-data reconstruction problem for tho 
space shuttla in used to demonstrate the applica- 
tion of the algorithm. Tho algorithm exhibits 
superior performance as compared with a standard 
implementation of the discrete Kalman filter. The 
algorithm has potential for solving many typos of 
nonatationary estimation problems for which the 
standard implementation of the Kalman filter is 
unsuited. 
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Figure 2, Topodatic axis coordinate definitions. 


Figure 1, Air-data estimation concept . 
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Figuro 6, Tina history comparisons of StS- 3 air data 
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crete Kalman filter. 
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